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ABSTRACT 

Aims. Most of current models of TeV blazars emission assume a Synchrotron Self-Compton mechanism where relativistic particles emit both 
synchrotron radiation and Inverse Compton photons. For sake of simplicity, these models usually consider only steady state emission. The 
spectral features are thus only related to the shape of the particle distribution, and do not depend on the timing of observations. 
Methods. In this letter, we study the effect of, firstly, the lag between the beginning of the injection of the fresh particles and the trigger of 
the observation, and secondly, of a finite injection duration. We illustrate these effects considering an analytical time-dependent model of the 
synchrotron emission by a monoenergetic distribution of leptons. 

Results. We point out that the spectral shape can be in fact very dependent on observational conditions if the particle injection term is 
time-dependent, particularly taking into account the effect of the time averaging procedure on the final shape of the SED. Consequences on the 
acceleration process are also discussed. 
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1. Introduction 

TeV blazars belong to the radio-loud active galactic nuclei 
(AGN) and are characterized by non-thermal continuum spec- 
trum, optical polarization, flat radio spectrum and strong vari- 
ability in all frequency bands. The broad band of the SED 
which extends from radio up to the TeV range, consists typ- 
ically in two bumps. The low energy one — peaking in the 
0. 1 - 1 00 ke V domain — is commonly attributed to synchrotron 
emission of ultra-relativistic leptons. In the Synchrotron-Self- 
Compton process (SSC) the second component is attributed to 
the inverse Compton (IC) scattering of synchrotron photon field 
by the same population of leptons. Emission of TeV photons 
provides evidence for the presence of particles at roughly the 
same energy if we suppose that they r adiate in the IC Klein- 
Nishina regime (iBlumenfhal & G ould 1970). Therefore, it re- 
quires an extremely efficient particle acceleration mechanism 
at work in the close environment of supermassive black hole. 

Extended obser vation s of TeV blazars like Mrk 501 (z s = 
0.034, lOuinn et"afl 1 19961) show that synchrotron component 
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can be very variable in the UV/X-ray range. The simplest ap- 
proach is the well-known homogeneous or one-zone modeling. 
It considers the emission of a relativistic moving blob of plasma 
filled by a tangled magnetic field, where ultra-relativistic parti- 
cle are injected and can cool freely (via both synchrotron and 
inverse Compton scattering process). Generally, in these mod- 
els, the acceleration mechanism is not consistently taken into 
account and the resulting energy distribution function (EDF) of 
injected particles is prescribed. If we supposed that the accel- 
eration zone could be dissociated from the radiative one, then 
we must distinguish two different kinds of EDF : the injected 
one and the cooled one. In the usual one-zone model, the first 
is directly the result of the acceleration mechanism. The sec- 
ond ensues from the radiative cooling of particles (and pos- 
sible re-acceleration mechanism). It is mathematically speak- 
ing the solution of the standard kinetic equation which the in- 
jected EDF is the main source term. The exact shape of the SED 
could be only understood in the framework of time dependent 
modeling, depend ing on the detail of the co oling and injection 
processes jKirk et alJl998llChiaberge & Ghiselfinil 19991) . The 
majority of SSC models are steady-state: therefore the instan- 
taneous SED is equal to the observed one. In this case, the 



2 



Ludovic Sauge and Gilles Henri: Time dependent homogeneous TeV blazars SED 



authors use broken power-law EDF as the solution of the ki- 
netic equation (characterized by the indices at and a u ). Then 
the subsequent synchrotron spectra have also a broken power- 
law shape 1 T v oc y~^ fj < with indic es 2 0f. u = (o7.» - l)/2 ^1/ 3 
over suitable frequency ranges jBlumenthal & Gould Il970). 
Spectral index variations can arise from the variation of ae, u 
i.e. by varying the physical con ditions of the emission zone 
(see e.g. lKino & Takaharal2004l) . 

We have recently investigated another solution, still in the 
framework of homogeneous modeli ng, but with a t ime depen- 
dent particle injection term JSauge & HenrilEoQ4l , hereafter 
paper I). The kinetic equation is solved numerically taking into 
account ( i) the possible in-situ pair reprocessing and ( ii) both 
synchrotron and IC cooling term. The source term of the equa- 
tion is chosen as a quasi-monenergetic distribution of particles 
(or "pileup"), injected during a finite time. Physically, the for- 
mation of such a distribution results from the combination of a 
stochastic hea ting via second order Fermi process an d radiative 
cooling (Henri & Pellet ierlll99lt lSchlickeiser 1985). Then the 
resulting cooled EDF is also partially a power-law, but with a 
constant index value a = 2 over a dynamical range depending 
on the details of the cooling process. In this case, the instanta- 
neous synchrotron spectrum is a power-law of index B - 1/2 
and no spectral index variation could be expected at this stage. 
However, such a variation can be obtained by averaging out the 
instantaneous SED over the time and considering that the ob- 
servation starting time does not necessarily coincide with the 
beginning of injection phase. 

In this letter we study in more detail this effect. First of 
all, we derive a basic model in Sect. 2 in the idealized case of 
monoenergetic injected EDF subject only to synchrotron cool- 
ing. Then we derive the global shape of the time average syn- 
chrotron part of the SED as a function of the observational 
parameters. Finally, we illustrate our approach on BeppoSAX 
archival data sets of the object Mrk 501 in Sect. 3, before con- 
cluding. 



2. The model 

In the following, all times are expressed in the plasma (blob) 
rest frame. We neglect the photon flight time across the source 
by considering evolution only less than R/c where R is the 
transverse size of the source (see Pap er I). This issue was p re- 
viously studied in detail by Chiaberge & Ghisellini ( 1999) in 
the case of the instantaneous emission. 



2.1. Energy Distribution Function of Particles 

For analytical purpose, we consider the simple case of the in- 
jection of a pure monoenergetic EDF of particle of energy 
7max,o during a finite time f inj 



<2inj(y; t) - Qq 6(y - 7max,o) n(f; 0, fi„j), 



(1) 



1 Only true in the case where particle energy dynamical range, i.e. 
the ratio y max /ymin is much larger than unity. 

2 At low frequency, this latter can not be steeper than the one pro- 
duced by a monoenergetic EDF, given the lower limit 1/3. 



and we consider that the main channel of particle cooling is the 
synchrotron process. We also neglect the in-situ pair creation 
process which can strongly affect the resulting cooled EDF in 
case of large yy-opacity (see Paper I). We introduced the gate 
function Yl(x; Xmi n , x max ) = ®(x - x m j n )0(x max - x) where is 
the usual Heaviside unit step function. We define the cooling 
time of a particle with a Lorentz factor y by the usual relation 
fcooi(y) = 1/ksyny, where k sya = cPx^/fmrriecB 1 ; It also repre- 
sents the time spent by an initial infinite energy particle to cool 
down to y. In the following, another characteristic time value 
related to the previous one will be t a \ t = f C ooi(7max,o)- Under 
these assumptions the solution of the standard kinetic equa- 
tion is analytic : n(y; t) = n y~ 2 n(y; y mm (f), y max (f)) where the 
power-law index 2 is typical of the synchrotron radiative cool- 
ing. Values of time-dependent bounds write 



, . ymax.O . . 

ymin(0 = "J — — ; ymaxW = 

1 + 1 1 1 cr it 



ymax.O 



1 +max(0,f-finj)/fcrit 



(2) 



We also define the dynamical range of distribution as ie the 
ratio r dyn (f) = y m ax(f)/ymin(0- 

2.2. Instantaneous synchrotron spectrum 

As said before, in the case where the previous dynamical 
range is sufficiently wide, the synchrotron spectrum is itself 
a power-law of index 1/2 in VT V . We recall that a particle of 
energy y emits preferentially synchrotron radiation roughly at 
the maximum frequency v s (y) = 0.29 x 3y 3 <y,(y)/47r where 
the relati vistic synchrotron gyro-freq uency writes u s {y) = 
qB/ym e c (Blu menthal & Gould ll970). Then we consider a 5- 
approximation of the synchrotron kernel function in the opti- 
cally thin emission limit, centred around v s . In this case, the 
synthetic synchrotron component SED can be approximated by 
VF v {t) = S v 1/2 n(v; (f)) over the frequency range 

given by 



x« = 



l(t) = 



"smax,0 



[l +max(0,f-f m j)/f cr it] 2 



(3) 



(1+f/fcrit) 2 ' 

where v smaXj0 = v smax (f = 0) = v s (y max , ). 

2.3. Time averaged synchrotron spectrum 

As previously said, the observed spectrum results from the time 
average of instantaneous spectra from f Q b s to f b s + Af b s . The 
time origin t — is related to the beginning of the injection of 
fresh particles. One writes, 



{vTv), = 



S v 



1/2 /"V„bs+A'obs 



A f ob: 

S v 112 

A?obs 



- dfn(v;v sl 

■s J/„ bs 

(Jl + Ii), 



n (0, 1 sin ax (0). 



(4) 
(5) 



where the quantities I\ and Ii arise from the splitting of the 
previous integral over [f bs»fcij] and [fmj,f bs + Af obs ]. After 
straightforward algebra, these values depend on the position of 
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f b s compared to f in j, and we distinguish two cases, 

• case A : f obs < finj 

I \ = max[0,f inj - max(f obs ,f crit (v))], 

1 2 = min[f obs + Af ob s,finj + fcrit(v)]-max[f in j,f crit (v)], 

• case B : f obs ^ fjnj 

Ji = 0, 

1 2 = min[f obs + Af obs ,f in j + f crit (v)]-max[f obs ,f cri t(v)]. 

where we define f crit (v) = f crit x [(v smax ,o/v) 1/2 - 1] representing 
the time needed by the lower bound of the synchrotron spec- 
trum to reach the frequency v. 

3. Results and Discussion 

3.1. Effects of parameters 

Using the previous relations, we show the influence of the f obs 
parameter on the time averaged spectra on Fig. [I] On theses 
figures, the time is normalized in unit of fj n j and frequency in 
unit of y S max,o- If fobs + Af obs > fi n j, the resulting SED consists of 
two parts (see left panel). The lower one exhibits a power-law 
tail of constant index 1 /2 and results from the emission of most 
energetic particle injected since the start of the observation (and 
not from the beginning of the injection) and which have not 
yet had the time to cool. Conversely, up to Vbreak the spectrum 
softens. In this case we obtain from Q v brea k = v smax (f ob s) or 
explicitly 

/ \ V S ma X ,0 

V brea k(fobs) = — : — rr. (O) 

(1 + fobs /fcrit) 2 

In the case f obs ^ fjnj ^ f obs + Af obs , the high energy part could 
be also described by a power-law with local index a, which 
implicitly depends on f obs . Combining previous calculations, 
we obtain 

log [<v!Tv)(Vs max ,o)/<v? r v)(v brea k)] 



log [v smax ,o/vb l eak J 

log [(l - f ob s/finj) (1 + fobs /fcrit)] 



21og(l + f obs /fcrit) 2 

where the upper bound is obtained for f obs 



(7) 



(8) 



0. Note that the 



spectrum is expected to be flat [a = 0) when f obs = f; n j - f a -it- 
This mechanism gives a plausible explanation of the broken 
power-law or curved s hape of TeV bla zars reported by many 
authors (see e.g. iFossati et al.ll2000albl) . Following our study, 
we deduce that the low energy synchrotron index should be al- 
ways be equal to 1/2 in VF V representation. Up to v brea k, the 
index value strongly depends both on the observational char- 
acteristics and on the details of the cooling process. Moreover, 
for fobs » fjnj - f cr it, this part of the SED strongly steepens and 
akin to a pure power-law. It is interesting to note that the max- 
imum curve naturally moves from v sm ax,o to v bre ak as the time 
increases and strongly depends on the value of observational 
parameters. Right panel shows the time evolution of spectral 
index calculated in three different logarithmic frequency bands, 
namely [-0.5,0], [-1,-0.5] and [-1,0]. 

Similarly, Fig. |2 shows the main effects of the increase of 
the Af obs parameter. In the case of f obs ^ fi„j ^ f obs + Af obs , the 





Normalized frenquency v/v syn (Hz) 



Fig. 1. Influence of the value of f obs parameter on the time av- 



0.5f 



inj- 



eraged SED (left panel) with Af obs = 4fi n j and f c , 
Vertical dashed lines represent location of the frequency break 
^break calculated for each previous curve (see eq. |6))- Resulting 
local spectral indices in normalized frequency range [-0.5,0] 
(upper solid curve), [-1, -0.5] (lower solid curve) and [-1,0] 
(dashed curve) are given on the right panel. 
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Fig. 2. Same as Fig. [0 for the influence of the Af obs parameter 
with f obs = f cl -i t = 0.5fj n j. Considering these parameter values, 
we easily verify that spectral index above v brea k = 0.25 v smaXj o 
is going to be flat since Af obs > f; n j - f obs = 0.5. 



spectral shape does not evolve. As expected, only the global 
flux level changes (decreasing linearly with Af obs ) just as the 
low-frequency cut-off (decreasing as maximum system time 
f obs + Af obs increases). 

3.2. Application to Markarian 501 

Mrk 501 is a privileged TeV blazar target. It has been stud- 
ied extensively during many multi-wavelength campaigns and 
have been observed in various spectral states. We now test our 
simple approach on state extracted from the BeppoSAX mis- 
sion public archive 3 . In the following we consider th ree differ- 
ent da ta sets, two from the 1997 April flaring period ( Pian et al. 
1998) and one from the 1999 June period (iTavecchio et alJ 
2001). Corresponding SED are given in Fig. [5] The two first 
SED have been studied in Paper I with the more complete ho- 
mogeneous SSC model and using high-energy data in order to 
constrain all physical quantities. 

All parameter values needed for the synchrotron part of the 
SED fits and suitable estimations are given in Tabled Again, 
the times are expressed in unit of f; n j . For each period, we rep- 
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Frequency v log[Hz] 

Fig. 3. Fits of Markarian 501 UV/X-ray synchrotron spectrum 
for three different periods. 



Table 1. Parameters used in our simulations. All times are ex- 
pressed in unit of f; n j 



data sample 


fcrit 




fobs 




Afobs 






min 


max 


est 




1997 April 7 


0.1095 


0.68 


0.92 


0.8 


0.946 


1997 April 16 


0.438 





0.5 


0.35 


4.2 


1999 June 


0.0512 


1.05 


1.15 


1.1 


4.5 



resent a continuous set of curves as the function of f b s ; also 
represented, our best estimation of f D b s . 

3.2.1. 1997 April period 

For the 1997 April 16 (medium state) and April 7 (high state) 
observations, we find of course values quite similar to Paper 
I. Fits are very good except for the extreme high energy part 
of the spectra. The exact form of this cut-off depends both on 
the exact profile synchrotron kernel and on the exact shape of 
the particle cooled EDF which are rather approximate in this 
work. For both periods, the beginning of the observation takes 
place before the end of the injection phase and therefore the 
synchrotron spectra displays a positive a value, as defined in 
equation |7}. 

First analysis done bv lPian et alJ dl998l) using an one-zone 
model with a power-law injection function show that the April 
16 spectrum is correctly reproduced for an injected EDF with 
an index 5=1 over a dynamical range less than 10. In the usual 
context of shock acceleration process, the value of the index s 
case, the index s directly depends on the physical c onditions 
at the shock (Blandford & Ostrikerl ll978l IProtheroe & Stanevl 
1999), 

face R+2 

5=1 + — = —— >1. (9) 

»esc K — I 

In the last expression, f acc and f eS c are respectively the acceler- 
ation and the escape characteristic time, and R the shock com- 



pression ratio. Then the value s = 1 is highly unlikely, because 
it requires f esc — > oo, which is not consistent : particles can not 
escape from the accelerator. Moreover, the narrowness of the 
dynamical range in this latter case clearly pleads for a monoen- 
ergetic injection function and consequently for a stochastic ac- 
celeration process. 

In classical shock acceleration model, the difference be- 
tween the 1997's high and medium state is explained by dif- 
ference in the shock physical condition — e.g. via the com- 
pression ratio R — leading to a different index value s. In the 
light of our scenario and as already noted in Paper I, the differ- 
ence in the shape of the spectra arises essentially from differ- 
ent observational conditions. The medium state corresponds to 
a previous injection observed in a later stage (with respect to 
beginning of the injection phase). This is corroborated by the 
flatness of the high energy part of the X-ray spectrum. 

3.2.2. 1999 June period 

During 1999 June period the synchrotron spectrum experienced 
a very steep state. In the context of the power-law injected EDF, 
this state is explained with a large value of the index a„ of the 
upper part of the EDF, roughly equals to 4.3 jTavecchio et alJ 
1200 lb . 

In the shock acceleration sc enario the value of ' a u is related 
to s via the relation a„ = s + 1 JKino & Takaharal2004h . Then 
in our precise case,it corresponds to an acceleration time quite 
larger than the escape one (f acc = 3.3 f esc or, in term of com- 
pression ratio R « 1.9). In our approach, this spectral shape 
is obtained considering an observation time ? b s larger than 

Anj — f crit ^ f inj ■ 

4. Conclusion 

We have exposed a simple time dependent mechanism in order 
to explain the spectral shape of the synchrotron spectrum of 
blazar considering the injection of a pure monoenergetic dis- 
tribution of ultra-relativistic particles over a finite time range. 
This latter EDF arises from stochastic acceleration mechanism 
(Second order Fermi process). Spectral index and SED shape 
variations can be explained by the value of the starting obser- 
vation time with respect to the beginning of the injection of 
fresh particles and from the variation of the acceleration con- 
dition even if the instantaneous synchrotron spectrum is uni- 
versal with a constant spectral index equals to 1/2. Because 
blazars are extremely variable sources, there are often observed 
in "Target-of-Opportunity" mode, where satellites observations 
are triggered by some other instruments after some delay lead- 
ing naturally to f b s + 0. This mechanism leads to complex time 
dependent behaviors of the spectral shape ; we show that we 
can reproduce various spectral shapes from the single to bro- 
ken power-law shape, with strong or soft cut-off. This model 
differs from the usual one in the context of the shock accelera- 
tion model, where the spectral variability arises from the phys- 
ical conditions at the shock. 
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